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f—( ' Abstract. Properties of neutrinos may be the origin of the matter- antimatter 

^ . asymmetry of the universe. In the seesaw model for neutrino masses this leads to 

important constraints on the properties of light and heavy neutrinos. In particular, 

an upper bound on the light neutrino masses of 0.1 eV can be derived. We review the 

present status of thermal leptogenesis with emphasis on the theoretical uncertainties 

^ ' and discuss some implications for lepton and quark mass hierarchies, CP violation and 

'nJ" ■ dark matter. We also comment on the 'leptogenesis conspiracy', the remarkable fact 

— »,^ I that neutrino masses may lie in the range where leptogenesis works best. 
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1. Introduction 

One of the great challenges of modern particle physics and cosmology is to explain the 
excess of matter over anti-matter observed in the universe. This baryon asymmetry is 
conveniently expressed as the ratio of baryon minus anti-baryon density to the photon 
density and has recently been measured to a high degree of accuracy by observations 
of the cosmic microwave background (CMB) pP combined with measurements of large 
scale structures of the universe j2] : 

In a complete cosmological model this baryon asymmetry has to be dynamically 
generated during the evolution of the universe in the hot and dense phase shortly after 
the big bang. This is possible if the particle interactions violate baryon number (B), 
charge conjugation (C) and the combined charge and parity conjugation (CP), and if 
the expansion of the universe leads to a deviation from thermal equilibrium |3] . 

All these ingredients are present in the Standard Model (SM) of particle 
interactions. However, the baryon asymmetry that can be generated in the SM falls 
far short of observations, i.e. an extended model of particle interactions has to be 
considered. The observation of neutrino masses also requires an extension of the SM 
and lepton number (L) violating interactions that are introduced in the seesaw model 
of neutrinos masses [H Ej can naturally give rise to the observed baryon asymmetry in 
the leptogenesis scenario [H] that is the topic of this article. 

The rather suprising fact that lepton number violating interactions can give rise to 
a baryon asymmetry is due to a deep connection between baryon and lepton number 
in the SM, as discussed in section |21 In section |3] we present the basic mechanism 
of leptogenesis and introduce some notations that are used in the following. The 
quantitative solution of the Boltzmann equations and the corresponding bounds on light 
and heavy neutrino masses are described in section 4. Here we also comment on the 
leptogenesis conspiracy, the remarkable fact that neutrino masses may lie in the range 
where leptogenesis works best. In section 5 the connection between low and high energy 
CP violation and implications of leptogenesis for the heavy neutrino mass spectrum are 
briefly discussed. Section 6 deals with the implications of leptogenesis for dark matter. 

2. Baryon and lepton number violation in the SM 

In the SM both baryon and lepton number are classically conserved, since they are 
protected by global Abelian symmetries. However, due to the chiral nature of weak 
interactions, these symmetries are anomalous and are violated at the quantum level 
[Zj. This is related to the non trivial vacuum structure of non- Abelian gauge theories, 
like the SM. Neglecting fermion masses, there are an infinite number of degenerate 
ground states whose vacuum field configurations have different topological charges, or 
Chern-Simons numbers Ncs |Hlini- In the electroweak sector of the SM a change in the 
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Figure 1. One of the 12-fermion processes which are in thermal equilibrium in the 
high-temperature phase of the standard model. 



topological charge, i.e. a transition from one vacuum to another one, corresponds to a 
change in baryon and lepton numbers, 

AB = AL = UfANcs , (2) 

where nj is the number of generations of quarks and leptons, i.e. n/ = 3 in the SM. 
Note that, although both baryon and lepton number are violated, the linear combination 
-B — L is still conserved at the quantum level. 

At low temperatures, when the electroweak symmetry is broken, the different 
vacua are separated by a potential barrier, whose height is determined by the vacuum 
expectation value (VEV) of the Higgs field, v = (0), i.e. the scale of electroweak 
symmetry breaking. Hence, processes changing the topological charge are tunneling 
processes whose rate is unobservably small, due to the smallness of the electroweak 
coupling constant. In the low temperature regime being probed at accelerator 
experiments, B and L are therefore conserved to a very good approximation, in accord 
with experimental observations (cf. [TU]). 

When the standard model particles form a heat bath of temperature T the situation 
changes. At high temperatures, T > Tew ~ 100 GeV, the Higgs VEV 'evaporates', 
leading to a restoration of the electroweak symmetry and the disappearance of the 
potential barriers separating the different vacua. B and L violating transitions are then 
no longer suppressed [TTj . 

The rate at which these processes occur is related to the free energy of field 
configurations which carry topological charge. In the electroweak part of the SM these 
so-called sphaleron processes lead to an effective interaction of all left-handed fermions 
^(cf. Fig.P), 

Ob+l = Yl {iLiququhi) , (3) 
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which indeed violates both baryon and lepton number by three units but conserves the 
combination B — L, in accord with Eq. Q. 

The sphaleron transition rate in the symmetric phase of the SM has been evaluated 
by combining an analytical resummation with numerical lattice techniques ^21- The 
result is that sphaleron processes are in thermal equilibrium for temperatures in the 
range 

100 GeV <T< 10^2 GeV . (4) 

These processes have a profound effect on the generation of the cosmological baryon 
asymmetry. Eq. Q suggests that any B + L asymmetry generated at temperatures 
T > Tew, will be washed out. However, since only left-handed fields couple to 
sphalerons, a non-zero value oi B + L can persist in the high-temperature, symmetric 
phase if there exists a non-vanishing B — L asymmetry. An analysis of the chemical 
potentials of all particle species in the high-temperature phase yields the following 
relation between the baryon asymmetry rjs and the corresponding L and B — L 
asymmetries tjl and tjb-l, respectively [Hj, 

Vb = cisph Vb-l = ^^j-r Vl ■ (5) 

Here a^ph is a number 0{1). In the SM with three generations and one Higgs doublet 

one has asph = 28/79. 

3. Leptogenesis 

3.1. A qualitative overview 

The deep connection between baryon and lepton number in the early universe has led 
to the realization that lepton number violating processes, whose presence is predicted 
by the seesaw model for light neutrino masses jHE], can be responsible for the observed 
cosmological baryon asymmetry. 

In the seesaw model the smallness of light neutrino masses is explained through the 
mixing of left-handed neutrinos with right-handed neutrinos i^r which are not present 
in the SM but are predicted in certain models of grand unification. The interactions of 
the SM are supplemented by the following Yukawa couplings of neutrinos, 

Cy = khuR(j) + u^MuR + h.c. , (6) 

where M is the Majorana mass matrix of the right-handed neutrinos, and the Yukawa 
couplings h yield the Dirac neutrino mass matrix mj:, = hv after spontaneous breaking 
of the electroweak symmetry. Since the Majorana mass matrix M is independent of 
electroweak symmetry breaking, one can have M ^ m^, which leads to the mass 
eigenstates 



z/ 



r-^ 



VL^yl = y\ N^VR^vl = N\ (7) 



m. = -mnijml . (8) 
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with heavy neutrino masses M and the hght neutrino masses 

1 
'M 

The heavy neutrinos are unstable and decay through their Yukawa couphngs into 
SM lepton and Higgs doublets. Due to their Majorana nature, the heavy neutrinos A^ 
can decay both into leptons and anti-leptons, i.e. lepton number is violated in these 
decays. In conjunction with B + L violating sphaleron transitions this leads to the 
required non-conservation of baryon number. Further, violation of C and CP comes 
about since the Yukawa couplings h are, in general, complex, thereby making possible 
the generation of a non-vanishing baryon asymmetry in these decays. 

A further complication is that the heavy neutrinos also mediate lepton number 
violating scattering processes which can erase any lepton asymmetry ^SlEI- However, 
the interaction rates for these processes are suppressed by the large mass of the heavy 
neutrinos if the temperature T is smaller than their mass. Hence, the lepton asymmetries 
produced in decays of the heavier neutrinos A^2,3 will be erased by lepton number 
violating scatterings mediated by the lightest of the heavy neutrinos, A'^i. Therefore, in 
the simplest case of hierarchical heavy neutrino masses. Mi -C M2,3, only decays of Ni 
can potentially explain the observed baryon asymmetry. 

The required deviation from thermal equilibrium is provided by the expansion of 
the universe. When the universe has cooled down to a temperature of order the heavy 
neutrino mass Mi, the equilibrium number density becomes exponentially suppressed. 
If the neutrinos are sufficiently weakly coupled they are not able to follow the rapid 
change of the equilibrium particle distribution once the temperature falls below their 
mass. Hence, the deviation from thermal equilibrium consists in a too large number 
density of heavy neutrinos compared to the equilibrium density JHl- Technically this 
requires the total decay width of Ni to be smaller than the expansion rate, the Hubble 
parameter H, at the time of decay, i.e. when T ~ Mi. This is the case if the effective 
neutrino mass, defined as 

ft. = ^^%^. (9) 

is smaller than the equilibrium neutrino mass 



m^ 



167r5/2 ^ 



V 



2 



where we have used Mp = 1.2 x 10^^ GeV and g^ = 106.75 as effective number of 
relativistic degrees of freedom in the plasma. The effective neutrino mass rhi is a 
measure of the strength of the coupling of Ni to the thermal bath. 

In order to see whether this mechanism of leptogenesis j^ can indeed explain the 
observed baryon asymmetry a careful numerical study is needed. As we shall see, 
successful leptogenesis is possible for rhi < m* as well as mi > m^,. A quantitative 
description of this non-equilibrium process is obtained by means of kinetic equations. 
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3.2. Boltzmann equations 

The evolution of particle number densities in the early universe is influenced not only by 
interactions but also by the expansion of the universe. It is convenient to scale out the 
expansion by considering the particle number Nx in some comoving volume element 
instead of the number density nx- For definiteness, we choose the comoving volume 
R*{t)^ which contains one photon at a time t* before the onset of leptogenesis, 

Nx{t) = nx{t) R,{t)\ (11) 

The final baryon asymmetry is expressed in terms of the baryon-to-photon ratio rjB, to 
be compared with the observed value rj^^^ (cf. Eq. (^). This is related to the B — L 
asymmetry in a comoving volume element by 

VB = ^iVi_„ (12) 

where / = 2387/86 is the dilution factor due to the production of photons from the 
onset of leptogenesis until recombination, assuming the standard isentropic expansion 
of the universe. 

Further, it is convenient to replace time t hj z = Mi/T, where Mi is the mass of 
the decaying neutrino. This is possible, since in a radiation dominated universe both 
variables are related by the expansion rate of the universe, the Hubble parameter H, 

ff = a = ,/i^Mi. (13) 

V 45 Mp ^2 ^ ^ 

In the simplest case of a hierarchical mass spectrum of right-handed neutrinos, 
Ml <^ M2,3, a numerical description of leptogenesis is provided by a set of two coupled 
differential equations |T7| ITHj. 



dNN^ 
dz 



(D + S) {Nx-,-K^^) , (14) 



^^^ = -6,D {Nx, - iV^^J - WNb-l , (15) 

where the terms on the right-hand side describe the effects of particle interactions. There 
are four classes of processes which contribute: decays, inverse decays, AL = 1 scatterings 
and AL = 2 processes mediated by heavy neutrinos. The term D accounts for decays 
and inverse decays, while the scattering term 5* represents the AL = 1 scatterings. 
Decays also yield the source term for the generation of the B — L asymmetry, the first 
term in Eq. (fTH|) . while all other processes contribute to the total washout term W which 
competes with the decay source term. Note that, in thermal equilibrium, Nx, = ^^, 
no B — L asymmetry can be generated, illustrating the need for a deviation from thermal 
equilibrium. The amount oi B — L asymmetry being produced by the source term is 
controlled by the CP asymmetry si in the decay of A'^i. 

In order to understand the dependence of the solutions on the neutrino parameters, 
it is crucial to note that the interactions terms D and S as well as the contribution from 
A^^i exchange to W are all proportional to the effective neutrino mass, 

D, S, ly. oc ^ , (16) 
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Figure 2. The expansion rate of tlie universe and the three interaction rates nomahzed 
to the expansion rate at z = 1 for a typical choice of parameters, Mi — lO^*' GeV, 
mi = 10^^ eV and m ~ 0.05 eV. 



whereas the strength of the remaining contribution to the washout term, AW, is 
determined by m^ = ml + rn^ + m|, the sum over the light neutrino masses squared. 



AW oc — ^ 



(17) 



If one assumes a vanishing initial B — L asymmetry before the onset of leptogenesis, 
i.e. at z <^ 1, the solution for Nb~l has the simple form 

Nb-l{z) = -^ ^1 f^{z- mi, M{fff) , (18) 

where we have introduced the efficiency factor k, |^ which does not depend on the CP 
asymmetry ei and parametrizes the effect of scattering and decay processes. It is given 
by the following integral expression: 



K,[Z 



{z) = ^ dz'D {N^,-N^^l)e-^^''^"^'^^'\ 



(19) 



It is normalized in such a way that its final value Kf = fi:(oo) approaches one in the 
limit of thermal initial abundance of heavy neutrinos, A^iv^(-2 ^ 1) = A^^^ = 3/4 and no 
washout, W = 0. 

As an example, the three interaction rates, (Td, Ts, T]y) = H z {D, S, W), and the 
Hubble parameter are shown in Fig.Elfor a typical choice of parameters. Mi = 10^° GeV, 
fhi = 10~^eV and fn = 0.05 eV. As can be seen from the figure, the out-of-equilibrium 
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Figure 3. The evolution of the A^i abundance and the B — L asymmetry for the same 
choice of parameters as in Fig. |21and ei = 10^^. 

condition fhi < 10"^ eV is fulfilled at 2; = 1, i.e. all interaction rates are smaller than 
the Hubble parameter. 

The corresponding evolution of the A^i abundance and the B — L asymmetry is 
shown in Fig.lHl starting ai z = 0.1 with a vanishing initial A^^i abundance. Although the 
neutrino production rates F^) and Fg are not in thermal equilibrium, they are still strong 
enough to produce a non- vanishing abundance of neutrinos at z < 1 and the equilibrium 
distribution is reached at z ~ 2. The required deviation from thermal equilibrium can 
clearly be seen as a small over- abundance of neutrinos for z > 2. However, the decay 
rate F^i also comes into thermal equilibrium at 2; ~ 2 leading to a rapid decay of the Ni 
abundance and the production of a non- vanishing B — L asjTiimetry. The change in sign 
in the asymmetry at 2 ~ 3 is due to the fact that the source term in eq. (|T5|) changes sign 
once the A^^i abundance becomes larger than the equilibrium abundance, i.e. neutrino 
production processes at 2; < 2 lead to a 'wrong sign' asymmetry that partially cancels 
against the asymmetry produced in the A^^i decays at z > 2. 



4. CMB constraints on neutrino masses 



We are now ready to discuss the implications of thermal leptogenesis for light and heavy 
neutrino masses. Combining Eqs. ()12|) and ()18p one obtains 

■r]B = —d El Kf (mi, MiTff) , (20) 
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where d = 3asph/(4/) ~ 10~^. Hence, determining the amount of baryon asymmetry 
produced in leptogenesis requires the calculation of both the final efficiency factor K{ 
and the CP asymmetry ei. A comparison with the observed value (cf. Eq. (P)) then 
allows to place stringent constraints on the involved seesaw parameters and, remarkably, 
on light and heavy neutrino masses. 

4.I. Final efficiency factor 

Starting from Eq. (fT^ and assuming a high initial temperature, Zi = Mi/Ti -C 1, 
the final efficiency factor can be calculated analytically pi^ I2()j. For values Mifff' -C 
10^^ GeV [m/(0.05eV)]^ the washout term /SW can be neglected and the final efficiency 
factor depends only on the effective neutrino mass rhi. 

The results for the efficiency factor are summarized in Fig. |3J Two different regimes 
can clearly be distinguished. In the weak washout regime, rhi <^ 10^'^ eV, the results 
strongly depend on the initial conditions and on theoretical assumptions, i.e. in that case 
predictions are strongly model dependent and affected by large theoretical uncertainties. 
On the other hand, in the strong washout regime, rhi ^ 10~^eV, the dependence on 
the initial conditions is practically negligible and the theoretical uncertainties are small 
such that the final asymmetry can be predicted within ~ 50%. 

In both cases very precise analytical approximations for the final efficiency factor 
can be obtained. It is instructive to start with a simplified picture where the scattering 
term S is neglected and only decays and inverse decays contribute. The decay parameter, 
^ ^ Td{z = oo) ^ mi 
H{z = 1) m* ' 

controls whether A^^i decays are in thermal equilibirium or not. Here To{z = oo) is 
the Ni decay width and m^, marks the boundary between the weak and strong washout 
regimes, as discussed in Section ITTl 

In the strong washout regime inverse decays rapidly thermalize the heavy neutrinos 
Ni and the washout due to decays and inverse decays is strong enough to destroy an 
initial asymmetry that may have been present before the onset of leptogenesis |21j . 
leading to a negligible dependence on initial conditions. Further, the integrand in 
Eq. (fT^ is peaked at a value zb ':^ I, implying that the final asymmetry is produced 
around a well defined temperature of baryogenesis Tb = Mi/zb ^ Mi, where the 
heavy neutrinos are fully non-relativistic. This contributes to reducing the theoretical 
uncertainties in the strong washout regime, since it has been shown that the Boltzmann 
equations employed here can be derived from a fully consistent quantum mechanical 
treatment in terms of Kadanoff-Baym equations in the non-relativistic limit j22] • 

In the strong washout regime the Ni abundance closely tracks the equilibrium 
abundance and a simple expression for the final efficiency factor can be obtained. 
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Assuming thermal initial A^i abundance, i.e. A^j^^ = 3/4, this expression also reproduces 
the correct asymptotical limit in the weak washout regime, K '^ 1. In Fig. |3] this 
analytical result is represented by the short-dashed line which has to be compared to 
the numerical results, given by the solid lines. 

In the weak washout regime the calculation is more involved if one starts from a 
vanishing initial Ni abundance. Indeed, in this case the Ni production by inverse decays 
leads to a negative contribution to the efficiency factor, corresponding to a 'wrong-sign' 
asymmetry, as discussed in Section 13.21 The final efficiency factor is thus the sum of a 
negative contribution, fi;j^(A'), and a positive one, k,^{K). A very good approximation 
for these contributions that interpolates between the weak and strong washout regimes 
is given by 

^^{K) = - 2 e-i^(^) (ei^(^) - l) , (23) 

4iK) = -^ (l - e-i-(^)^^W) . (24) 

Here N{K) ~ (9'n-/16)K is the maximal Ni number density being produced in the 
weak washout regime and N{K) interpolates between N{K) and the maximal number 
density A'^eq = 3/4 in the strong washout regime. For large K the negative contribution 
is suppressed, while the positive one asymptotically approaches Eq. (j221)- On the other 
hand, in the weak washout regime the positive and negative contributions cancel each 
other to leading order in K, i.e. the total efficiency factor is of order K^ [2S|, as shown 
in Fig. m where the short-dashed line again corresponds to the analytical solution and 
the solid one to the numerical integration of the Boltzmann equations. 

This cancellation of the leading order contributions to the final efficiency factor 
in the weak washout regime no longer occurs when the scattering term 5* is taken 
into account. Indeed, these scattering processes enhance the A^i production thereby 
giving rise to an additional positive contribution to the efficiency factor. On the other 
hand, these scatterings are CP-conserving, i.e. they do not contribute to the negative 
part of the efficiency factor, as long as the contribution of AL = 1 scatterings to 
the washout term is negligible which is always the case in the weak washout regime. 
In this way scatterings can greatly enhance the final efficiency factor in the weak 
washout regime |21]. The drawback is that the result is very sensitive to different 
approximations being used in the computation of the scattering rates and different 
results, ranging from the case where scatterings are negligible to a behaviour K{ oc K, 
have been obtained. Potentially important effects that have recently been discussed but 
are presently controversial include scattering processes involving gauge bosons [23 121] 
and thermal corrections to the decay and scattering rates ^Tf [22] • The range of different 
results is represented in Fig. |3]by the hatched region. An additional uncertainty in the 
weak washout regime comes about due to the dependence of the final results on the initial 
A^i abundance and a possible initial asymmetry created before the onset of leptogenesis. 

The situation is very different in the strong washout regime. The final efficiency 
factor is not sensitive to the neutrino production since a thermal neutrino distribution 
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Figure 4. Final efficiency factor when the washout term AW is neglected. 



is always reached at high temperatures, z <^ 1. However, AL = 1 scatterings also 
contribute to the washout term but their effect is small and thus the theoretical 
uncertainty arising from these scattering processes is not larger than about 50%. This 
uncertainty for values mi > m^, ~ 10"^ eV is again indicated by the hatched region 
in Fig. El A much more important source of uncertainties are spectator processes [ I28j . 
which can change the produced asymmetry by a factor of order one. 

In the strong washout regime the highest efficiency is reached when scatterings 
are neglected and only decays and inverse decays are taken into consideration, which 
approximately corresponds to the results obtained in j^EI- In this parameter range the 
final efficiency factor is given, within theoretical uncertainties, by the simple power law 



Kf = (2 ±1)10" 



0.01 eV 
mi 



l.litO.l 



(25) 



The only model independent information we have on mi is that it has to be larger than 
the smallest neutrino mass mi [221 • However, a situation where mi < mi <^ Trisoi is 
rather artificial within neutrino mass models and the leptogenesis predictions are then 
very model dependent in the weak washout regime. In typical neutrino mass models 
values of mi are usually in the mass range suggested by neutrino oscillations. It is 
remarkable that both the scale of solar neutrino oscillations, ttIsoI = A/Am-^^j ^ 8 
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10 eV, and the scale of atmospheric neutrino oscillations, m 
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are much larger than the equihbrium neutrino mass m^,. Hence, the parameter range 
suggested by neutrino oscillations, msoi ^ "2i < matm, hes entirely in the strong washout 
regime where theoretical uncertainties are small and the efficiency factor is still large 
enough to allow for successful leptogenesis. 

4-2. CP asymmetry 

The CP asymmetry ei is the second crucial ingredient needed to calculate the baryon 
asymmetry. To leading order in the Yukawa coupling /i, the CP asymmetry is determined 
by the interference between tree level and vertex plus self-energy one-loop diagrams 
JSOIEII and can be consistently extracted from the Icj) — >■ 14> scattering processes [32|- 

We will be interested in the maximal CP asymmetry for given neutrino masses. 
Assuming a hierarchy of the heavy neutrinos, M2^z ^ ^ii ^^^^ depends on Mi, mi, mi 
and m^ [21J. It can be expressed as 



£r'(Mi,mi,mi,m3) =£r'(Mi)/3(mi,mi,m3) , /3 < 1 . (26) 

it is given by 



The maximal asymmetry, i.e. /3 = 1, is reached for mi = and, with 7713 = \/m1^^ + ?7if , 



3 Mimatm ,^-6 f ^h \ f "iatm 



,„„„„ - ^ ^^^^^ - 10- [j^^,^) (^^j . (27) 

An improved bound is obtained in the limit nii/rhi -^ 0, where one obtains [SI], 

^(„0 . !!!l^:!!l . (28) 

Also important is the case of quasi-degenerate neutrinos, m3 ~ mi. In this region one 

finds isniEn], 



/3mi,mO^^ 1 /l_^. 29 

For all neutrino mass models with moderately hierarchical heavy neutrinos, such that 
Eq. (|26|) applies, the observed baryon asymmetry yields constraints on the three neutrino 
parameters mi, mi and Mi. In the following two sections we follow mostly the discussion 
in Ref. [Hj. 

4-3. Lower hounds on heavy neutrino masses and reheating temperature 

The maximal baryon asymmetry r]]^'^^ [Mi, rhi, mi) is the asymmetry corresponding to 
ef"^^ (Ml, mi, mi). The CMB bound then amounts to the requirement 

r7r^(Mi,mi,mi)>r/g^^^. (30) 

This represents an interesting constraint on the space of the three parameters Mi, mi 
and mi. We have seen that the absolute maximum of the CP asymmetry is obtained for 
mi = 0. For rui 7^ the function (3 suppresses the CP asymmetry [33]. Furthermore, the 
AL = 2 washout term becomes stronger when the absolute neutrino mass scale increases 
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(cf. (J17j) ). Therefore, the maximal baryon asymmetry r7^'^^(Mi,mi,?7ii) is maximized 
for nil = 0, and in this case the allowed region in the space of the parameters Mi and 
•fhi is maximal [THj . 

In this way one finds an important lower bound on the value of Mi [331 EI; by 
inserting the expression (pTjl into the CMB constraint (jHn|l (cf. fj2n|l ). 

1 16 7r v^ T^g*^^ 



d 3 matm l^i 



Vb^'^ \ f0.05eV\ _i 



Ml > M^"^ = 

~ 6.4 X 10«GeV ( -^ ] ( ^^^^^^^^ ) t^V' . (31) 

For thermal initial abundance, and in the limit rhi/m,, — > 0, one has by definition /tj = 1, 
and therefore 

Ml > (6.6 ± 0.8) X 10^ GeV > 4 x 10^ GeV . (32) 

Here the last inequality is the 3 a bound, with the experimental value of Eq. (P) for 
rf^^^ , and matm = (0.051 ± 0.004) eV. In the case of a dynamically generated Ni 
abundance the maximal efficiency factor is K{ ~ 0.18, which yields the more stringent 
bound 

Ml > (3.6 ± 0.4) X 10^ GeV > 2 x 10^ GeV . (33) 

The most interesting case corresponds to the range rrisoi ^ "^i ^ ""^atm, for which 
the power law Eq. ()25|) for Kf can be used. Using the central value of kj and neglecting 
the theoretical uncertainty, one obtains 

Ml > (3.3 ± 0.4) X 10^° GeV (^F^) ' (^4) 

which yields the 3 a bound 

Ml > (1.5 - 10) X 10^° GeV . (35) 

The lower bound on Mi is particularly interesting since it can be translated into a lower 
bound on the initial temperature Tj which, within inflationary models, corresponds to 
the reheating temperature. 

So far we have assumed that the temperature T[ is larger than Mi. If one relaxes 
this assumption, the flnal efficiency factor is in general reduced. For small rhi, however, 
the threshold value for T^, below which the reduction is appreciable, is given by Mi 
itself. Below this temperature the A^^i abundance is either Boltzmann suppressed, for 
thermal initial abundance, or the A^^i production is considerably suppressed, for zero 
initial abundance. Therefore, for values fhi < 10^'^ eV, the bounds (|3^ and ^^ apply 
also to the reheating temperature Ti. 

In the more interesting case of strong washout, about 90% of the flnal baryon 
asymmetry is produced in a temperature interval zb — 2 < z < zb + 2. Hence 
the reheating temperature can be about {zb — 2) times lower than Mi without any 
appreciable change in the predicted flnal asymmetry [20] . In the interesting range 
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^soi ^ "^1 ^ "^atm one has 2;^ — 6 — 8, and therefore the bound Eq. (j^ gets relaxed 
by a factor 4 — 6, such that 

Ti > (4 X 10^ - 2 X 10^°) GeV . (36) 

Compared to the small mi range, mi < 10^'^ eV, the lower bound on the reheating 
temperature is slightly more restrictive in the favoured range mgoi ^ ^1 ^ "^atm due to 
the loss in efficiency. 

4-4- Upper bound on light neutrino masses 

For large values of the absolute neutrino mass scale the AW washout term cannot be 
neglected. The final efficiency factor can be calculated using the approximation that 
AW starts to be effective for z > zb, where the asymmetry generation from decays has 
already terminated. This works very well in the strong washout regime. Since mi > mi, 
this does not introduce any restriction if m-i > m* ^ 10~^ eV. One then has 

_ _ " ( ,yi "] (~\^ 

Ki{mi, Mifn^) ~ Kf(mi)e "s Uo^jgcvAcv; ^ (37-) 

where u ~ 0.186 and /t/(mi) is the efficiency factor calculated in the regime of small 
neutrino masses, neglecting the AW term. Because of the assumption of strong washout 
one can use the simple power law Eq. ()25j) . At certain peak values of Mi and mi the 
maximal baryon asymmetry rj]^'^^ reaches the absolute maximum 

with X — 1-6 eV'^ and ^ oc m^^.^^^. This finally yields the leptogenesis bound on the 
neutrino masses, mj < 0.1 eV [21^. Deviations of the quantity C, from unity account 
for a change of the input parameters as well as various corrections, such as a possible 
enhancement of the CP asymmetry or supersymmetry. 

A more precise calculation has to take into account the dependence of neutrino 
masses on the renormalization scale. The running of the atmospheric neutrino mass 
scale matm from the Fermi scale fi = mz to the high scale n ^^ Tb makes the bound 
less restrictive. On the other hand the bound on m^ obtained at high energies has to be 
evolved down to low energies. This second effect is dominant and thus taking the scale 
dependence into account makes the neutrino mass bound more restrictive jHZj- The 
smallest effect is obtained for a Higgs mass Mh — 150 GeV, which leads to a ~ 20% 
more restrictive bound P^. Taking into account this effect one then obtains the Sa 
bound 



m 



, <0.12eVe'/'. (39) 



Note, that the bound m^ < 0.1 eV conservatively accounts for the theoretical 
uncertainties including spectator processes [2H1 which make the bound more restrictive 
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by about 0.02 eV |. The strong suppression of the baryon asymmetry with increasing 
neutrino mass scale, tjb oc 1/m^, which is reflected in the dependence $,^^^, makes the 
bound rather stable. This is different from the lower bounds on Mi and Tj which relax 

as 1/^ mi. 

It is important to keep in mind that the neutrino mass bound can be evaded, 
with some effort. A measurement of the neutrino mass scale above the leptogenesis 
bound would require significant modifications of the minimal leptogenesis scenario that 
we described. The possibilities include quasi-degenerate heavy neutrinos J2I1 ES] , non- 
thermal leptogenesis scenarios jHHl IHH HHj or a non-minimal seesaw mechanism as in 
theories with Higgs triplets [HI US US! • 

4-5. Leptogenesis conspiracy 

The upper bound on the neutrino masses ()39|) arises when the information on matm 
from neutrino mixing experiments is employed. The value of rriatm sets the scale for 
the transition from a hierarchical, with mi <^ ^^atm? to a quasi-degenerate neutrino 
mass spectrum with mi > "matm- The joint action of the CP asymmetry suppression 
for mi > matm and rhi ~ m,i (cf. (1221) ), together with the washout from the AW 
term, place a limit to the level of degeneracy of the light neutrino masses and, using the 
measured value of the atmospheric neutrino mass scale, an upper bound on the absolute 
neutrino mass scale. 

We now want to study how this upper bound gets relaxed if the experimental 
measurement of the atmospheric neutrino mass scale is ignored. Since the maximal 
asymmetry is obtained for hierarchical neutrinos, i.e. /? = 1, we have to use the simple 
bound (|?r|). This also implies that the upper bound on the absolute neutrino mass scale 
will coincide with an upper bound on the atmospheric neutrino mass scale itself, since 
rn ~ rriatm — ^^3- The maximal baryon asymmetry is then approximately given by 

r]T''{Mi,mi,ms) = der^'iMi) Ki{mi,Miml) . (40) 

Using Eq. pTjl it is easy to see that the maximum of the asymmetry is realized for 

MiC:,2zBlO''GeY ^^l^i^V , (41) 

and is then given by 

r7^^"(mi, ma) ~ 0.7 x 10"^ Zij(mi) Kf (mi) ( 2:2^ ) . (42) 

\ rn3 J 

I In Ref. [2SI the upper bound 0.15 eV has been obtained, which is 0.03 eV weaker than the bound H39|l . 
About 0.02 eV of this difference is due to the different treatment of radiative corrections which depend 
on the top and Higgs masses. The remaining 0.01 eV reflects differences in the treatment of thermal 
corrections. This is included in the theoretical uncertainty of the efficiency factor Kf(mi) (cf. Eq. H25|) '). 
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In the case of zero initial A^i abundance the maximum is obtained for mi ~ 
2 X 10~^eV, where k{ ~ 0.2 and 2;^ ~ 2 §, implying 

vT'-ims) - 3 X lO-^' (^^^^) • (43) 

Note that the peak lies in the strong washout regime where results do not depend on 
the initial conditions. 

For rris ~ m^.tm — 0.05 eV, Eq. (jlH)) is in good agreement with the numerical results 
of Ref. ^H]- If we do not make use of the experimental information on rriatm and just 
require that the peak asymmetry is larger than the observed value given by Eq. (^, 
then we obtain the upper bound 

ms < 250 eV . (44) 

Together with Eq. ()4H) this implies the lower bound for heavy neutrinos masses, 

Ml > 2 X 10^ GeV . (45) 

This exercise shows that without the experimental knowledge of matm the bound on light 
neutrino masses would have been much looser. However, it also demonstrates, even more 
remarkably, that the neutrino oscillation data, together with the laboratory bounds on 
light neutrino masses, represents a highly non trivial test of thermal leptogenesis. 

As discussed in the previous section, the leptogenesis bound of 0.1 eV appears 
to have a theoretical uncertainty of about 0.03 eV. For comparison, ten years ago 
it was believed, based on the same equations, that Majorana masses m^ ~ 10 keV 
and Ml ~ 1 TeV were compatible with thermal leptogenesis [Ej. During the past 
years theory and experiment contributed about equally, on a logarithmic scale, to the 
progress. On the theoretical side the better understanding of the Boltzmann equations 
was important whereas the measurement of the atmospheric neutrino mass scale was 
the crucial experimental ingredient. 

Actually, the successful matching of thermal leptogenesis predictions with 
experimental data is even more intriguing. From Eq. ()42|) we can derive an upper bound 
on fhi by using the strong washout behaviour Kf ~ 2/{zbK) (cf. (J22|l ) and imposing 



the CMB bound, which yields, with matm = A/m| — 



m 



l5 



mi<20eV('^:^^^) , (46) 

'^atm 



again in very good agreement with the numerical results [18]. From this expression 
one reads off that only for matm ^ leV it is possible to have mi ~ matm- If the 
experiments had found matm ^ 1 eV, thermal leptogenesis would have worked only 
for models with mi <^ frii^tm, requiring a large amount of fine tuning, as already 
pointed out. Furthermore, the favoured strong washout regime implies matm > fn*- 
Hence, leptogenesis favours for the atmospheric mass scale the range 10~^ eV — 1 eV, 
in remarkable agreement with experimental data. Thanks to this 'conspiracy' the 

§ This can be inferred from Fig. 4 and Fig. 5 of Ref. j2()| . 
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seesaw mechanism, for the same values of the involved parameters, explains equally 
well both the neutrino masses and the observed baryon asymmetry, with a remarkable 
independence on the assumptions about the inflationary stage or, more generally, the 
cosmological stages that precede leptogenesis. 

This conspiracy would become even more impressive if the lightest neutrino mass 
should turn out to be larger than m^. This would imply, in a completely model 
independent way, rhi > m*, with thermal leptogenesis in the strong washout regime. 
Together with the upper bound mi < 0.1 eV, this selects the optimal leptogenesis window 
10^^ eV < nil ^0.1 eV for the absolute neutrino mass scale. 

5. Flavour aspects 

Leptogenesis is closely related to other processes involving neutrinos and charged 
leptons. The upper bound on the light neutrino masses and the lower bound on the heavy 
neutrino masses have already been discussed in the previous section. Very interesting 
are also the connection with neutrino mixing and with lepton flavour changing processes 
in supersymmetric theories. During the past years these subjects have been studied in 
great detail by many groups (cf. jlU EHl HHl HIl ESI)- In the following we shall discuss 
two important examples, the possible connection between CP violation in neutrino 
oscillations and leptogenesis, i.e. at low and high energies, and the relation between 
leptogenesis and the heavy neutrino mass spectrum. 

In the standard model with right-handed neutrinos the masses and mixings of 
leptons are described by three complex matrices, 

'Cm = Wmoi^R + eLmiCR + v^Mur + h.c. . (47) 

For the Majorana mass matrix M one expects M ^ mo-, which leads to three light and 
three heavy Majorana mass eigenstates, z/ ~ z/j^ + i/£ = i/'^ and A^ ~ z/^j + z/^ = A^'^. 

In the following we will work in a basis where M is diagonal and real, with 
Ml < M2 < M3, which is appropriate for leptogenesis. The light neutrino mass 
matrix ttIj, and the charged lepton mass matrix mi are then diagonalized by the unitary 
transformations, 

f mi \ 

V^''^^m^V^''> = - m2 = mf^^ , (48) 

\ ms y 

/ me \ 

V^'^^meV^'^ = m^ = mf'^ . (49) 

I m^ J 

The leptonic (MNS) mixing matrix 

U = v(^)tvM (50) 

describes the couplings of mass eigenstates in the charged current, 

4l = -jfLruuLW- +... , (51) 
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which leads to neutrino oscillations. 

It is well known that, in general, CP violation in neutrino oscillations and in 
leptogenesis are unrelated [IH]. The reason is that the CP asymmetries Ei in heavy 
neutrino decays only depend on m/j^mD. Hence, changing rriD to Ktud, where K 
is a general unitary matrix, leaves Ei invariant whereas the leptonic mixing matrix 
U = ^(^)t^('^) ig changed to U = V^^^^ KV'^'^\ and therefore arbitrary. Still, the question 
remains whether in some physically well motivated cases a connection between CP 
violation at low and high energies exists. 

CP violating observables are most conveniently described by weak basis invariants 
which are inert under a unitary transformation, / — > Kl, of the lepton doublet 
/ = {vli^l}- For neutrino oscillations the appropriate variable is the commutator 
between the hermitian matrices Hy = mym^ and Hi = miml |5lH l5T] , 

TT[H,,Hif OC A,2lA,32A,3iAi2lAi32Ai3iJi, (52) 

A^2i = ml-ml,..., Ai2i = ml - ml , . . . , (53) 

where J/ is the leptonic Jarlskog invariant j52j . 

Ji = Im[UnU22Ut2U;i] oc sin(5 , (54) 

which is proportional to the CP violating phase 6 of the mixing matrix U. In a basis 
where the charged lepton matrix mi is diagonal and real one has 

TT[H,,Hif OC Ai2iAi,2Aisilm[H,,2H,23H,si] ; (55) 

correspondingly, for diagonal and real m,,y one has 

Tr[iJ„ Hi]' ex A,2iA,32A,3i lm[H 1,2^23^^1] . (56) 

The CP asymmetry in Ni decays can be conveniently expressed in terms of the 
weak basis invariant |22] 

El OC lm.[m'^m,„m,*j-,]ii . (57) 

Comparing Eqs. (J56|) and (J57j) the independence of CP violation at low and high energies 
is obvious. In a basis where m,^ is diagonal and real, CP violation in neutrino oscillations 
is entirely determined by the phases m; whereas the lepton asymmetry only depends on 
the phases of m^. 

An interesting example is the case of only two heavy Majorana neutrinos, A^^i and 
N2, which corresponds to the limit where A^s decouples jSHj- For a given texture of m^ 
one then obtains (cf. (jHSj), (jFfj) ). 

El OC lm[H^i2H^23Hu3i] OC sin (5 . (58) 

Hence both, the CP violation at low and at high energies, are determined by the Dirac 
phase 6. A further low energy quantity is the Majorana phase entering neutrinoless 
double beta decay. In some models with maximal atmospheric neutrino mixing this 
phase coincides with the leptogenesis phase jHl]. All CP phases can also be related in 
models of spontaneous CP violation [5111 133]. 
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Another important question is the connection between leptogenesis and the heavy 
neutrino mass spectrum. In the lepton mass eigenstate basis the neutrino mass matrix 
can be reconstructed from data, i.e. the leptonic mixing matrix and the neutrino masses, 

m^ = -U^^'^mf^'^U^"'^'^ . (59) 

The heavy neutrino mass matrix is then determined by 111^ and tti^j, 

M-^ = -m^^m^{m^)-^ . (60) 

Using the SO(IO) mass relation mn = m„, where rriu denotes the up-type quark mass 
matrix, and assuming that in the lepton mass eigenstate basis rriu is also diagonal and 
real, rriu = diag(mu, mc,mt), one obtains the heavy neutrino masses in terms of rriy, 
iriui '"^c and irit jHEI- Neglecting the hierarchy among the small neutrino masses one 
estimates that the hierarchy among the heavy neutrinos is very large. 



^ ^ :::!i ^ iq-^" . (6i) 

M3 \mt) ^ ' 

With M3 ~ m^/yA^^ ~ 10^^ GeV, this implies Mi ~ 10^ GeV. A detailed study 
yields masses in the range Mi ~ 10^ . . . 10^ GeV. As we saw in the previous section, 
such small masses are incompatible with conventional thermal leptogenesis. A possible 
way out are quasi-degenerate heavy neutrino masses jSHl- Alternatively, the most naive 
SO(IO) mass relations are not correct. 

Parameters consistent with thermal leptogenesis have recently been obtained in a 
six-dimensional SO(IO) GUT model, compactified on an orbifold jHZI. Due to mixings 
with a heavy lepton doublet and a heavy right-handed down-type quark one finds for 
the mass matrices in a particular flavour basis the relations, 

M oc m„ , m'^, ~ m'^ ~ m'^ . (62) 

Here M and m„ are diagonal 3x3 matrices, and m^, m'^ and m'^ are 4x4 matrices, due 
to the mixing with the heavy states. Integrating them out, the neutrino mass matrix 
can be expressed approximately in terms of the quark mass matrices, 

niy oc nid mT . (63) 

For the light neutrino masses this leads to the estimate 

!^~f!^V!^~0.1. (64) 

A detailed calculation |^ confirms this result. Furthermore, one finds the mixing angle 
©13 ~ {mjn},) / {nitms) ~ 0.1, and e ~ 10"^, rhi ~ 10~^ eV, Mi ~ {■mu/mt)M3 ~ 
10^'^ GeV, consistent with thermal leptogenesis. 
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6. Cosmological aspects 

It is well known that the temperature required by thermal leptogenesis, Tb ^ 10^ GeV, 
is potentially in conflict with the thermal production of gravitinos in the early universe 
|58| l59] . Late gravitino decays after nucleosynthesis significantly alter the successful 
BBN predictions. 

The production of gravitinos is dominated by QCD processes, and the gravitino 
number density increases linearly with the reheating temperature after inflation, 

where g^ is the QCD gauge coupling. Correspondingly, the BBN upper bound on 
the allowed gravitino energy density, P3/2 = 'm3/2n3/2, implies an upper bound on 
the reheating temperature. Detailed studies jHOJ IHI] lead to the stringent bounds 
Tr < 10^ — 10^ GeV for gravitino masses in the range 777,3/2 = 0.1 — 1 TeV. Hence, 
unstable gravitinos are in conflict with thermal leptogenesis, unless the gravitino mass 
is very large, r?73/2 ^ 50 TeV. 

Non-thermal leptogenesis models are still compatible with the above bounds on 
the reheating temperature. For instance, in some supersymmetric models the scalar 
partner A^i of the heavy neutrino A^^i might be the infiaton [B2] and its decays could then 
generate the baryon asymmetry. In this case the leptogenesis temperature can be below 
10^ GeV jHHl I2E], which would be consistent with the above gravitino bounds on Tr. 
However, a recent analysis of the BBN constraints [H^ || with particular attention to the 
hadronic decay modes of the gravitino yields the much stronger bound Tr -C 10^ GeV 
for gravitino masses 7773/2 = 0.1 — 1 TeV. Hence, unless the gravitino is extremely heavy, 
also non-thermal leptogenesis appears to be inconsistent with unstable gravitinos. 

Already in connection with thermal leptogenesis, it has therefore been suggested 
that the gravitino may be the lightest superparticle (LSP) and stable jHS]- In this case, 
gravitino production is enhanced [HE] . 

n^ M^ \m3/2j 

where m-g is the gluino mass. Consistency with BBN and the observed amount of dark 
matter then yields an upper bound on the gravitino mass and a lower bound on the 
mass of the next-to-lightest superparticle (NSP). In [02] the case of a higgsino NSP 
was analyzed, which is now disfavoured due to the improved BBN bounds on hadronic 
NSP decays jHl]. Still viable is the case where a scalar lepton is the NSP jHZI. For 
Tb = 3 X 10^ GeV the upper bound on the gluino mass is m-g < 1.3 TeV if a scalar tau 
is the NSP; for a scalar neutrino NSP one finds m-g < 1.8 TeV j67j . 

The relic density of gravitinos is determined by thermal production at the reheating 
temperature Tr after inflation and also by NSP decays after their freeze-out temperature. 
It is an interesting possibility that the latter process dominates. This would be the case 

II Note that the analysis strongly depends on the assumed prhiiordial ^Li abundance. 
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Figure 5. Relic gravitino density for different values of reheating temperature and 
gravitino mass, ^jrf' = 1. m-g = 1 TeV, which implies m.3/2 < 0.1 TeV for a stable 
gravitino. For Tr > T*, ^-^jiY? is independent of Tr and 7713/2- 



for low reheating temperatures, incompatible with leptogenesis. One then obtains a 
prediction for the amount of gravitino dark matter which is independent of the reheating 
temperature [HH] ■ If gravitinos are the only component of dark matter, the superparticles 

1 TeV and vtir^ > 0.5 TeV 



have to be rather heavy. For t^ as NSP one finds 7/13/2 = 0.2- 



Recently, it has been pointed out that the thermal production of gravitinos is 
significantly changed in theories where the gauge coupling depends on the expectation 
value of a scalar field ^ fTU]. For instance, in the case of gaugino mediation one has, 

1 



1 0T 
M 



9l 



9'{<Pt) 



(67) 



where Qq is the zero-temperature gauge coupling, M is a mass scale between the 
unification scale and the Planck mass, and 0t is the deviation of the field from its 
zero-temperature value at temperature T. At temperatures above a critical temperature. 



T* ~ -m 



3/2 



M. 



m~g 



1/2 



(68) 



the gauge couplings decreases, and the gravitino production is frozen. Remarkably, the 
relic gravitino density is essentially determined just by the gluino mass [70] . 

m~g \3/2 /^^/^ 



^3/2/1^ - (0.05 - 0.2) X (- 



TeV/ 



T]^ 



(69) 



where the factor ^/rf ^ depends on the mechanism of supersymmetry breaking. For 
gaugino and gravity mediation one has ^/r/^ = 0{1). The observed amount of dark 
matter is then obtained for a gluino mass m-g ~ 1 TeV, which will be tested at the LHC 

^ Here the parameter ^^' ^ denotes the ratio of dilaton mass and gravitino mass. 
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As Fig. 5 illustrates, the temperature where thermal leptogenesis takes place is likely 
to be larger than the critical temperature T*. For a gluino mass rrig ~ 1 TeV one then 
obtains automatically the observed amount of cold dark matter. 

Finally, thermal leptogenesis can also be consistent with gravitino dark matter if 
the gravitino is very light, 7713/2 — 0.1 — 10 MeV [7T| 17^ . which is realized in gauge 
mediation models. Alternatively, the gravitino can also be very heavy, ms/2 ~ 100 TeV 
|73| I74j ■ as in anomaly mediation. 

In summary, leptogenesis strongly constrains the nature of dark matter. For 
very heavy unstable gravitinos, 7713/2 ~ 100 TeV, where ordinary WIMPs can be the 
dark matter, the superparticle mass spectrum is strongly restricted. Alternatively, the 
gravitino has to be the LSP with a mass below ~ 0.1 TeV; the observed value of QcDAih"^ 
can then be naturally explained. Hence, the identification of the invisible dark matter, 
which will hopefully take place at colliders in the coming years, will also shed light on 
baryogenesis and therefore on the origin of the visible matter. 
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